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It is well known that boundary conditions on quantum fields produce divergences in the renor- 
malized energy-momentum tensor near the boundaries. Although irrelevant for the computation 
of Casimir forces between different bodies, the self-energy couples to gravity, and the divergences 
may, in principle, generate large gravitational effects. We present an analysis of the problem in 
the context of quantum field theory in curved spaces. Our model consists of a quantum scalar 
field coupled to a classical field that, in a certain limit, imposes Dirichlet boundary conditions on 
the quantum field. We show that the model is renormalizable and that the divergences in the 
renormalized energy-momentum tensor disappear for sufficiently smooth interfaces. 

PACS numbers: 03.70.+k;11.10.Gh 

I. INTRODUCTION 

The vacuum energy produces measurable forces between neutral bodies Q]. But even in the presence of a single 
body, the modes of the electromagnetic field are disturbed and produce a self-energy. Although irrelevant for the 
calculation of forces between different objects, this self-energy is in principle observable through its coupling to gravity. 

In the case of perfect conductors, the vacuum energy density, or more generally the energy-momentum tensor, 
diverges near the boundaries, as noted for the first time by Deutsch and Candelas a long time ago These 
divergences are not the usual ones in quantum field theory, because they are present in the already renormalized 
energy-momentum tensor. The origin of the divergences is the unphysical assumption of perfect conductivity for all 
modes of the electromagnetic field since, on physical grounds, one expects any material to become transparent at high 
energies. More generally, divergences in the renormalized energy-momentum-tensor are present even for non-perfect 
conductors, as long as there is a sharp boundary between two media with different electromagnetic properties. In this 
case the origin of the divergences is that, for modes with extremely small wavelengths, it is unphysical to assume a 
sharp boundary, since the transition region between media becomes larger than the wavelength of the high frequency 
modes. 

The problem of the divergences in the self-energies has been considered in flat spacetime in a series of works by 
Graham et al [3]. It was pointed out that, if the boundary conditions are replaced by an interaction with a second 
field, the model is renormalizable and the self-energies are not divergent. In those works, the authors considered a toy 
model consisting of a vacuum scalar field cj> coupled to a background (classical field) tr, with an interaction of the form 
<f) 2 a. The sharp limit corresponds to consider a discontinuous a, and the "perfect conductor" limit to take a — > oo 
at a particular point, imposing Dirichlet boundary conditions on the vacuum field at that point. When these limits 
are taken, the self-energy depends on the ultraviolet cutoff, i.e. on the specific interaction between the high energy 
modes of the quantum field and the microscopic degrees of freedom of the body. Recently, Milton [4] and Bouas et al 
[5] considered a similar problem, computing the energy density for a scalar vacuum field in some particular potentials 
("soft walls"), showing that the energy density near a potential barrier is finite for sufficiently smooth potentials. For 
a discussion of some aspects of the coupling of the Casimir energy to gravity see Ref.[B]. 

A complete analysis of the Casimir self-energies and their eventual gravitational implications must be performed in 
the context of quantum field theory in curved spacetimes. Following Ref.[3], one could replace the boundary conditions 
by interactions with a classical field, and show not only that the matter sector of the theory is renormalizable, but 



* fdmazzi@cab.cnea. gov. ar 
' jnery@ic.sunysb.edu 
t alcsatz@umd.edu 



2 



also that the usual divergences in the energy-momentum tensor can be absorbed in the coupling constants of the 
gravitational sector, resulting in finite and well defined semiclassical Einstein equations when the classical field is 
sufficiently smooth. The divergences associated to the unphysical limits will reappear when considering discontinuous 
classical fields. The aim of the present paper is to provide such analysis. We will consider a scalar vacuum field in 
curved spacetimes, coupled to a classical field a that models the "mirror". We will show that, for smooth a, the 
model is renormalizablc using the standard renormalization procedure for quantum fields in curved spacetimes. We 
will provide some examples to illustrate the appearance of divergences for sharp interfaces. 

It is worth remarking that similar divergences in the renormalized energy-momentum tensor do appear even for free 
fields in curved spacetimes, if one assumes that the spacetime metric is not sufficiently smooth. A well known example, 
in cosmology, is the divergence in the energy density that appears when one considers an abrupt transition between 
two different epochs (e.g. De Sitter space to radiation dominated universe [7]), that involves a discontinuity in the 
curvature tensor. A less known example is the divergence of the vacuum polarization around spherically symmetric 
objects, at a sharp interface between the object and the vacuum [S]. 

The paper is organized as follows. In the next section we introduce the model. In Section 3 we prove that the 
divergences in the energy-momentum tensor can be absorbed into the bare constants of the theory, yielding finite 
semiclassical Einstein equations. In Section 4 we present some explicit calculations for weakly coupled mirrors, showing 
that as long as the potential that models the mirror is smooth enough, there are no divergences in the renormalized 
energy- momentum tensor and in the renormalized coincidence limit of the two point function, (4> 2 ). On the other 
hand, discontinuities in the potential or in its first two derivatives produce infinite answers. Section 5 contains a 
discussion of the main results. We will work with natural units K = c = 1 and metric signature H . 



II. THE MODEL 



We consider a quantum vacuum field <f> interacting with a background classical field a on a curved spacetime. The 
action of the complete system is 
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Here R^ = i? a Mal /, k = 87rG, A is the cosmological constant, and e^, i=l ,2,3 are dimensionless parameters. The 
terms quadratic in the curvature are needed in order to renormalize the theory, as is the self-interaction term for a. 
Note that the classical field a provides a position-dependent mass for the field cf>, and therefore the propagation of </> 
will be suppressed in regions where the mass is very high. In this sense a models a "mirror" . For example, a thin 
mirror located at x — xq is described by the interaction 



Ai<7 4> = Xi5(x — xq)4> 



(4) 



and the perfect conductor limit corresponds to Ai — > oo. 
The classical field equations are 
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where 



2 fiSgrav 

7^ 5 9^ 



-G 

K 



(8) 



and 



2 SSmat _ j,( a } 



(9) 



The tensors -ff^J come from the variation of the terms quadratic in the curvature contained in the gravitational action. 
The classical energy-momentum tensor for the a field is 



T$ = (1 - 2Z 2 )<T ;l ><T;p + (26 - l)g^; p a' P - 2^(7(7.^ 



6V 



(10) 



while Tfit' is the energy-momentum tensor for a free field with variable mass m 2 + \o 2 , that is, 

T$ = (l - 2£ 1 )0 ;/ > ; „ + (26 - Dg^-A'" - 2&# !|u ,+ 



(11) 



We consider now the semiclassical version of the theory, in which the field <f> becomes a quantum field while a and 
<7 M „ are treated classically. The Heisenberg equation for the quantum operator associated to <fi is given by the classical 
Eq.fpj). The evolution equations for the classical backgrounds are obtained by taking the mean value of the classical 
Eqs^6]) and Q: 



(12) 



(□ + m2+£ 2J R+^(0 2 )) ( 7+^7 3 = O (13) 
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We have therefore established a well defined model to study the divergences of the self-energy in Casimir calculations. 
The quantities (0 2 ) and {T^f) are formally divergent. The divergences must be absorbed into the bare constants of 
the theory. The renormalized version of (T^t ) contains information about the Casimir effect (force between different 
objects), as well as the self-energy. The gravitational effects produced by the vacuum energy can in principle be 
computed by considering it as a source of the semiclassical Einstein equations. 



III. RENORMALIZ ABILITY 



The theory of interacting fields in curved spacetimes can be renormalized using a precise covariant procedure [3]. 
In the present model it is necessary to introduce minor modifications to take into account that not only the metric 
but also one of the interacting fields is treated classically. 

It will be particularly useful to adapt the renormalization method described in Ref.[TU]. It is shown there that to 
analyze the renormalizability of Ay 4 theory in curved spaces at the level of the equations of motion, one can split 
the field as tp — ipo + <p, where ipo is the mean value of the field and (p is the quantum operator that describes the 
fluctuations around the mean value. To one loop order, ip satisfies a free field equation with a variable mass. The 
situation in our model is very similar, since the quantum field <j> can be thought of a free field with variable mass. 
Therefore, we can follow closely Ref . |10| . 
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We define the renormalized quantities 



•ul2 



X^/ii/ /ren V^ujy / an la, 



(14) 
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where (T^)adA and (</> 2 ) a d2 are constructed using the Schwinger DeWitt expansion up to fourth and second adiabatic 
order respectively [IT]. The divergences present in these quantities are to be absorbed into the bare constants of the 
theory. 

The usual Schwinger DeWitt expansion for the propagator of a massive field reads [T2] 



G ( s ] D {x,x') = -2lmG s F D (x,x') = 



= -21m 
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where G F is the Feynman propagator, A(x,x / ) is the van-Vleck determinant, a(x,x') is one half of the geodesic 
distance between x and x' and n is the number of spacetime dimensions, that acts as a regulator. The functions 
Q,j(x,x') are defined by a set of recursive equations that follows from imposing the equation for the propagator. 
When the field has a variable mass, this expansion can be generalized to [TO] 



G { sl(x,x') = -2ImG F (x,x') 
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ds 
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where f(x,x') = |[M 2 (;c) + Af 2 (x')] with Af 2 (x) = to 2 + ^a 2 . 

The quantities to be subtracted to cancel the divergences of (</> 2 ) and (T^) are [TP] 
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We emphasize that we are not using a point-splitting regularization but a dimensional rcgularization. We will substi- 
tute (T^t ) and (4> 2 ) in the semiclassical equations by (T^) ren + (T^P) ad4: and (</> 2 ) ren + {4> 2 ) a d2i respectively, and 



absorb the infinities contained in (2) 
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and 
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into the bare constants appearing in those equations. 



A. Renormalization of the equation for the classical field 



The coincidence limit of the two point function for a field with variable mass is given by [10] 
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where /i is an arbitrary constant with dimensions of mass. Inserting this expression into the semiclassical equation 
for <7 and writing the bare constants in terms of renormalized ones 



we obtain 
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Therefore, the divergences can be absorbed with the counterterms 
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where Am|, A£ 2 and AA 2 are finite contributions (they vanish in the minimal prescription scheme). Note that in the 
conformal case (mi = and £i = 1/6), only a counterterm for the self-coupling of a is needed. 



B. Renormalization of the semiclassical Einstein equation 



Let us now consider the renormalization of the gravitational sector of the theory. As n — > 4, the divergent part of 
(T$)adA is given by [B] 
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This expression contains both geometric divergences and divergences dependent on the classical field a (and its 
derivatives). The former should be absorbed into a redefinition of the gravitational constants appearing in the left 
hand side of the semiclassical Einstein equations, and the latter into a redefinition of the constants associated to the 
a field, which appear in T\ 



to 



Inserting Eq.(23) into Eq.(12| one can show that the divergences dependent on a can be absorbed using the same 
counterterms given in Eq.(22). This is a non trivial check of our calculations, and a necessary condition for the 
renormalizability of the theory. On the other hand, the geometric divergences can be absorbed into the gravitational 
constants, by choosing the following counterterns 
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(25a) 
(25b) 



This completes the proof of the renormalizability of the model, and is one of the main results of this paper. We have 
shown that, if the presence of a mirror is modeled by the interaction of the vacuum field with a classical background 
field, the divergences in the vacuum expectation value (T^f , ) can be absorbed into the bare constants appearing in the 



G 



semiclassical Einstein equations. Not only the gravitational constants are renormalized but also the bare constants 
associated to the Lagrangian of the classical background field. The renormalizability is valid for one or more mirrors 
of arbitrary shape, as long as they can be described by a smooth function a. Note that the divergences in the 
energy- momentum tensor are considerably simpler in the conformal case m\ = and £1 = 1/6. 

As a final remark, we stress that the quantities (0 2 ) ren and {T^t )ren will be finite for a sufficiently smooth 
background field a and spacetime metric g^ v . We will illustrate this in the next section. 



IV. EXPLICIT EVALUATIONS IN THE WEAK FIELD APPROXIMATION 

In this section we present simple expressions for evaluating explicitly ((f) 2 (x)) and (T$)(x)) for a given background 
potential c(x), within a weak field approximation. For simplicity we will work in Minkowski space, and set m\ = 
so that the quantum field is massless. The sense in which we define the weak field approximation is that we require 
M 4 <C D(M 2 ), or in other words X 2 a 4 (x) <C AiD(cr 2 (a;)). This can be achieved by having a weak coupling Ai, or a 
weak and/or rapidly varying field a(x). The calculation will stay at the lowest order in Ai. 

The procedure we follow is similar to the one used in |5] . It is based in solving the equation for the Feynman Green 
function, 

n+^-o 2 (x)^G F (x,x') = -8{x-x'), (26) 
to first order in Ai, and using the resulting G^p [x, x') to compute (4> 2 (x))^ and (T^P (x)} W in the coincidence limit. 



Expanding Eq.(26l by writing G F = G { p } + G { p in a Ai expansion and discarding the second-order term, we obtain 



□ G { P { Xl x') = -^-a 2 (x) Gf (x, x') . (27) 



This is solved explicitly using a momentum-space representation for Gp : 

\ r r r —ik(x — x) —ik'(x—x') 



Gp (x, x 1 ) = ^J**J ^ J — 2 . (28) 

After switching variables to p = k — k' and q = k + k' , the fc-integrals can be done by passing to rt-dimensional space 
and employing standard dimensional regularization techniques, with the result in the coincidence limit being: 

G ™<^> - -I^SVW - ^ I ^/dW*-V(*)ln (- J) . ,2 9 , 

Here [i is an arbitrary mass scale introduced in the regularization procedure. The divergence in the first term, which 
is purely local, can be absorbed in the coefficient £2 of the action, as can be seen by comparison with the second line 



of Eq.(22|. Hence in accordance with Eq.(17) we have 



(0 2 (x)) ro „ = ^ J d 4 p J d 4 xc^-^a 2 (x) In (-^j . (30) 

For example, in the particular case in which the background potential depends only on the spatial coordinate z, we 
shall have: 
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(31) 



This expression can be used to compute (4> 2 (z)) ren for the case of one or several parallel flat mirrors, which are "almost 
transparent" in the sense that the expression is valid to first order in the coupling Ai. 

We will prove now that according to this model (4> 2 (z)) ren is everywhere finite if the background potential <r 2 (z) 
is an intcgrable C 2 function. Under this assumption, it follows as a corollary of the Riemann-Lebesgue lemma that 
the Fourier transform of <r 2 (z) (which we note by a 2 (p z )) falls off faster than p~ 2 . Hence g(p z ) = ^n\p z \a 2 (p z ) is 
integrable, from which it follows that its Fourier transform is in turn well-defined and finite at all points. Since 



according to (31 1 it is precisely the Fourier transform of g(p z ) that gives the nonlocal part of (</> 2 (z)) ren , and the local 
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part is the term involving ln(/i 2 )tr 2 (z), it follows that both the local and the nonlocal parts of (4> 2 (z)) rcn are finite for 
all z if cr 2 (z) is an integrable C 2 function. 

We have therefore proved that the divergences in (<fi 2 ) are removed if the background field modeling the mirrors 
is sufficiently well-behaved. We are not able to give a general proof using weaker assumptions than C 2 continuity. 
However, we will see that in simple concrete examples the assumption of C° continuity is sufficient to obtain finite 
results. 

Let us first consider a discontinuous mirror of width 2z , specified by the background potential: 



a 2 (z) = 



1 



2Lz 



e(zo-M): 



(32) 



where Q(z) is the unit step function. The parameter L has length dim ension and is introduced so that <r 2 (z) is 
normalized by J o 2 {z) dz = L^ 1 . Introducing this expression into Eq.(31 1 we obtain as a result: 
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[sgf> + z ) ln(|z + z Q \ + 7) - sg(z - z Q ) ln(\z - z \ + 7)] 



(33) 



where we omit the local contribution coming from the ln(/i 2 ) term. Notice that the result diverges logarithmically at 
the boundary of the mirror, where <J 2 (z) is discontinuous. 
Consider next a continuous mirror given by 



a 2 (z) 
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z < — (z Q + d), 

— (z + d) < z < -Zq, 

- z < z < z , 
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if 
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1 if 



(34) 



wherein the discontinuities in Eq. ( 32 ) are resolved with linear interpolations of width d. For this mirror the vacuum 
polarization is computed to be: 



<^(*)>rc 
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(35) 



where we have omitted again a local contribution proportional to a 2 . This result is everywhere finite; it can be 
compared to that of Eq.(33) in the plot in Figure 1. 

We turn now to the study of (T$ (x)) to first order in Ai. From the expression (JTTJ) it follows that (T^ u ) (we drop 
the label from now on) is obtained from the coincidence limit of Gf(x,x') as: 



-Im 



(1 - 2Z 1 )d IM d , v G F (x,x') + (2^ - Vltu rrd p d>,G F (x,x>) 



Ai 



2 6 d^Gpix, x') + 2^ W pa d p d a G F {x, x 1 ) + -^ v ^G F (x, x') 



(36) 



The coincidence limit of the last term is ~ X\((j> 2 ), and since {(j) 2 ) is 0(\\) we may discard this term in our first-order 
approximation. The coincidence limit of the remaining term can be performed in an analogous way to that of {4> 2 ), 
by passing to n-dimensional space, using the momentum space representation (28) for the first-order Gf(x,x'), and 
using standard dimensional regularization techniques. The divergences appearing are proportional to o 2 niv +r) )Jil/ \3<j 2 , 
and can be absorbed into a redefinition of £2 as we did before. After the renormalization is performed and we return 
to n — 4 we are left with 



Ai 



5127T 6 



d 4 p 



d 4 xe ip ^-^a 2 (x)hx (-^j 



(37) 



which is analogous to (30) for (4> 2 (x)) rcn . For the case in which a 2 depends only on z we get by analogy with (31 ): 

Ai 



(T^(z)) Te 
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(38) 




FIG. 1: A comparison of {(j) 2 (z)) ren for a discontinuous potential (red) and a continuous one (blue), with z — 4 and 
d = 1. Note that (4> 2 {z)) r en diverges at the points where the background field a is discontinuous, and that the 

divergences disappear when a is continuous. 



Comparing with our arguments showing finiteness for (<fi 2 (z)) ren it is clear that the components of (T Miy (z)) rcn will 

not diverge as long as h(p z ) = p 2 ln\p z \a 2 (p z ) is integrable. To ensure this it is sufficient to require o 2 {z) to be an 
integrable C 4 function, since then its Fourier transform falls off faster than p~ 4 . However, for the simple particular 
example of a mirror similar to ( 34 ) but with polynomials with higher degree of continuity interpolating between zq 
and zo + d, we have checked that C 2 continuity of <J 2 (z) is enough to render the result finite. 

We stress at this point that the results obtained in this section have their gravitational counterparts: as shown in 
Rcf.[8], when one computes (4> 2 (z)} ren and (T Ml/ (z)) rcn for a free quantum field in a curved background, both quantities 
diverge at the points where the background is not sufficiently smooth. The examples described in Ref.[B] refer to 
the vacuum polarization around spherically symmetric objects, and divergences show up for sharp interfaces, where 
the matter density is discontinuous. The divergences disappear when the matter density and its first two derivatives 
are continuous across the interface. Similarly, in a cosmological context, divergences in the renormalized energy- 
momentum tensor are removed when one replaces an abrupt transition of the scale factor by a smoother transition 
in which the scale factor and its two first derivatives are continuous. The concrete example of the transition between 
the inflationary period and radiation domination is discussed in detail in Ref. [7j. 

High vacuum energy densities in the presence of boundaries may potentially produce large gravitational effects. In 
our model, the usual divergences are rendered finite when the sharp boundary is replaced by a sufficiently smooth 
background field. However, the quantum vacuum energy does attain high values if the background field varies over 
short distances. We c an make a crude estimation of the maximum value of the energy density using dimensional 
analysis. From Eq.(38), it is not difficult to see that if <J 2 (z) is varying over distances of order d then we should expect 
the maximum value of {T^ v {z)) lcn to be of order \i<7 2 /d 2 , where <7q is a representative value attained by <J 2 (z). The 
quantity AiCTq has units of (length) -2 , and determines the properties of the mirror. Therefore it will be identified 
with the square of the plasma frequency lo p . Restoring h and c factors, the mass density associated to the vacuum 
fluctuations is of order huj 2 /c 3 d 2 . If we take d to be in the range of the I0~ 10 meters, representing a smearing of 
the conductor's sharp boundary over atomic length scales, and uj p to be in the range of 10 15 Hz, a typical value of 
the plasma frequency, then we see that firstly, the first-order approximation assumed in this section is validated, and 
secondly, the peak values of the quantum energy densities at the boundary of the conductor are of order 10 _13 g/cm 3 , 
too small to be gravitationally detectable. 
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V. CONCLUSIONS 

In this paper we analyzed the coupling of the vacuum self-energy to gravity. We considered a model in which the 
presence of the bodies that disturb the modes of the quantum fields is described by the interaction with a background 
classical field. We have shown that the divergences in the energy-momentum tensor of the quantum fields are consistent 
with the semiclassical Einstein equations, that is, they can be absorbed into the bare constants of the theory. As 
expected, the divergences in the renormalized energy-momentum tensor noticed in previous works only appear when 
considering unphysical limits of perfect conductivity and/or sharp interfaces, that in our model would correspond 
to take a non-smooth background field. We have shown that modeling the mirrors in the Casimir effect with a 
background potential a 2 removes the divergences attached to the boundaries, as long as a 2 is a sufficiently smooth 
integrable function (C 2 continuity is sufficient for (<fi 2 ) rcn and C continuity for (T Ml/ ) rcn , though weaker assumptions 
are enough in simple examples). The proof is carried to first order in the coupling Ai- These results are in tune with 
those found in Refs.jUJS], where finiteness of the results is shown with exact calculations for some particular examples 
of <t 2 , and are analogous to those in Refs.[71[S] for the vacuum polarization of free fields in curved backgrounds, when 
the metric and its first two derivatives are continuous. We expect that similar conditions for finiteness can be obtained 
beyond the perturbative approximation. 

We have considered a toy model for a vacuum scalar field. We have presented some specific examples to illustrate 
the appearance of divergences for sharp interfaces. The examples suggest that the gravitational effects of the self- 
energy near mirrors are extremely small for realistic values of the relevant parameters. We expect similar results to 
be valid for more realistic models involving the electromagnetic field. 
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